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Abstract. We present results of a VLBA multi-frequency study of the 
pc-scale twin jet in NGC 1052. We observed this object at epoch 1998.99 
with the VLBA at 5, 8.4, 22 and 43 GHz both in total and linearly polar- 
ized intensity. The spectral analysis confirms the necessity of a free-free 
absorbing medium, obscuring the innermost part of both jets. At 5 GHz 
we found a compact linearly polarized emission region at the base of the 
eastern jet with a degree of polarization of 1.5%. At higher frequencies 
there is no evidence for polarization in our data. A core shift analysis 
constrains the position of the central engine to ~ 0.03 pc. The shift rates 
of the apparent core position with frequency confirm the strong influence 
of free- free absorption in conjunction with steep pressure gradients at the 
bases of both jets. 



Introduction: NGC 1052 is a nearby elliptical galaxy that harbors an active 
galactic nucleus in its center. Its radio structure ranges from ~ 20" down to 
the smallest resolvable scales below 1 milliarcsecond. Its bright compact radio 
core and its small distance of D = 22.6 Mpc makes NGC 1052 a prime target for 
VLBI studies of the environment of the central engine believed to be a super- 
massive black hole of about 1O 8±1 M . 

On VLBI scales NGC 1052 shows a jet /counter- jet structure oriented at a posi- 
tion angle of about 65°. Kinematical studies at 15 GHz Q show that the two 
jets are oriented close to the plane of the sky. The distribution of water maser 
emission in NGC 1052 on pc-scales was imaged with VLBI by Q. The maser 
emission is concentrated in two groups of unresolved maser spots, lying along 
the western jet. Strong free- free absorption at the corresponding region of the 
source was found by || (see also The presence of atomic gas at the cen- 
ter of NGC 1052 is evident from HI absorption observations Q and the X-ray 
spectrum of the nucleus ||. 



Imaging and model fitting: In Fig. [T] we present VLBA images of NGC 1052 
in total intensity at 5, 8.4, 22 and 43 GHz obtained from our observations on 
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Figure 1. VLBA images of NGC 1052 at 5, 8.4, 22 and 43 GHz. The 
map parameters are given in Table |. 

December 28th 1998 (related parameters in Table |]). The images have been 
produced following standard procedures. In the images we define four regions, 
A (eastern jet), B, C and D (western jet). The alignment of the images at the 
four frequencies was performed through a pairwise identification of Gaussian 
model fit components between adjacent frequencies. We derived spectra (shown 
in Fig. ^) of the different regions of jet and counter-jet separately, by adding up 
the flux densities of the corresponding model components. 

Both the eastern and the western jet have a steep spectrum at the highest 
frequencies. At low to intermediate frequencies the eastern jet remains steep 
between 22 and 8.4 GHz and becomes flat between 8.4 and 5 GHz whilst the 
western jet exhibits a turnover of its spectrum below 22 GHz. The innermost 
part of the western jet (region B) has a highly inverted spectrum. No emission is 
detectable in this region at 5 GHz, making the spectral index a (S oc u +a ) larger 
than 3. This value exceeds the theoretical maximal value of 2.5 for synchrotron 
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Figure 2. Spectra of the four source regions (A, B, C and D in Fig. |l|), 
obtained from the model fitting of the visibility data. 

self-absorption (see e.g. @) marking a possible region of external absorption (in 
an obscuring torus?) at the base of the western jet. 

Polarimetry: We carried out a polarimetric imaging analysis following stan- 
dard procedures (see J7j). We detected linearly polarized emission above a 3a 
level only at 5 GHz. At this frequency, we found a compact region of linearly po- 
larized emission at the base of the eastern jet, with a flux density of about 3 mJy 
beam" 1 . This corresponds to a degree of polarization of 1.5 %. The electric field 
vectors are parallel to the jet axis, i.e. the corresponding emission region is dom- 
inated by a magnetic field perpendicular to the jet axis (assuming no Faraday 
rotation). Our spectral measurements indicate that the emission from this part 
of the eastern jet is optically thick at 5 GHz. This implies that the polarized 
flux should increase at higher frequencies, a trend that we do not observe in our 



Table 1. Map parameters for Fig. [I]. 
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3.30 X 1.31 , - 


-3.74 


0.660 b 


2.41 


0.26 


0.99 


x (2,4,. ..,256,512) 
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1.98 x 0.81 , - 


-3.87 


0.538 b 


2.39 


0.25 


0.81 


x (2,4,. ..,256,512) 


22 


1.3 


0.86 x 0.32 , - 


-7.63 


0.339 c 


1.51 


1.20 


1.02 


x (2,4,. ..,128,256) 


43 


0.7 


0.45 x 0.16 , - 


-7.93 


0.126 c 


0.67 


0.67 


1.89 


x (2,4,...,16,32,64) 



a Total flux density recovered in the map model; b Corresponds to the A component; c Corresponds to 
the B component. 
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Figure 3. Left panel: 5 GHz image of the total intensity (contours) 
and the linearly polarized emission (segments representing the EVPA, 
proportional to the flux density). Right panel: profile (dashed line in 
the left panel) of the total and linearly polarized intensity at 5 GHz 
taken along P.A.= 69°. 

data. This can be interpreted as an effect of optical depth: if the jet is linearly 
polarized only in the outer layers, no linearly polarized emission is expected at 
high frequencies where the radiation originates in deeper layers of the jet, closer 
to the geometrical jet axis. 

Core-shift analysis: The symmetry between the jet and the counter-jet con- 
strains the position of the central engine in NGC 1052. The location where a 
jet becomes visible at a given frequency is usually referred to as the "core". 

The core is located at the distance to the central engine r c , where the optical 

i 

depth r has fallen to be ~ 1. This distance is given by: r c ~ v k r , where 
k r = ((3 — 2a)m + 2n — 2)/(5 — 2a), with m and n the power indices of the 
magnetic field and the particle density: B ~ r~ m , N ~ r~ n ||. Measuring r c 
at two frequencies allows one to determine k r in the corresponding region of the 
jet. For a freely expanding jet in equipartition ||, k r =l. The value of k r is 
larger in regions with steep pressure gradients and may reach 2.5, for reasonable 
values of m,n (<4, see [||). If external absorption determines the apparent core 
position, comparable density gradients of the external medium can alter k r to 
values above 2.5. 

The values of k r deduced depend crucially on the absolute values of r c on the two 
sides and therefore on the assumed position of the central engine. Four scenarios 
have been tested with different reference points (see Fig. |3]). Table || gives the 
derived values of k r for each scenario. The area between the model components 
A 15 and B 2b is the most likely location of the central engine and the center 
between both components is a natural choice for its exact position (scenario 1). 
Shifting the reference point eastwards (scenarios 2 and 3) alters the values of k r 
into unphysically large regimes (requiring density gradients oc r -10 and higher). 
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Figure 4. Image with the model fitting results of the 43 GHz data. 
The circles represent the Gaussian components. The innermost jet 
components are labeled as A 14, A 15 and B 2b. The putative loca- 
tions of the central engine are indicated (dashed arrows) for 4 different 
scenarios (see discussion in the text). 



Assuming the true center of activity to be located more westwards (closer to 
B2b, scenario 4), the values of k r derived are still acceptable. We show the 
core positions of both jets at the different frequencies for the first case (scenario 
1) in Fig. [|. The eastern jet has rather high values of k r below 22 GHz, although 
still in agreement with steep pressure gradients in the jet environment. Above 
22 GHz k r is 3.9=L0.8, which is a good indicator for free-free absorption affecting 
the jet opacity. The western jet has values of k r as high as 6.8±2.7 between 22 
and 43 GHz, suggesting a large contribution from free-free absorption. 
Our results from the core-shift analysis support the picture of a free— free ab- 
sorbing torus covering mainly the inner part of the western jet and also a smaller 
fraction of the eastern jet. The true center of activity in NGC 1052 can be deter- 
mined to lie between the model components A 15 and B 2b, with an uncertainty 
of only ~ 0.03 pc. 
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Table 2. k r values for the four different putative centers of activity. 
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Figure 5. Core positions in the two jets at the different frequencies 
for scenario 1 (see Fig. The corresponding values of k r for the 
core shift between adjacent frequencies are labeled. Smaller core shifts 
correspond to higher values of k r . Table [2] provides the values of k r for 
the other three scenarios. 
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